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Abstract. We analyze the dynamics of the entanglement in two independent non- 
Markovian channels. In particular, we focus on the entanglement dynamics as a 
function of the initial states and the channel parameters like the temperature and 
the ratio r between loq the characteristic frequency of the quantum system of interest, 
and u> c the cut-off frequency of Ohmic reservoir. We give a stationary analysis of the 
concurrence and find that the dynamic of non-markovian entanglement concurrence 
C p (t) at temperature fc^T = is different from the fc^T > case. We find that 
"entanglement sudden death" (ESD) depends on the initial state when ksT = 0, 
otherwise the concurrence always disappear at finite time when fc^T > 0, which 
means that ESD must happen. The main result of this paper is that the non- 
Markovian entanglement dynamic is fundamentally different from the Markovian 
one. In the Markovian channel, entanglement decays exponentially and vanishes only 
asymptotically, but in the non-Markovian channel the concurrence C p (t) oscillates, 
especially in the high temperature case. Then an open-loop controller adjusted by the 
temperature is proposed to control the entanglement and prolong the ESD time. 
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1. Introduction 

Entanglement is a remarkable feature of quantum mechanics, and its investigation is 
both of practical and theoretical significance. It is viewed as a basic resource for quantum 
information processing (QIP) p], like realizing high-speed quantum computation [2] and 
high-security quantum communication j3]. It is also a basic issue in understanding the 
nature of nonlocality in quantum mechanics [H [5j [6]. However, a quantum system 
used in quantum information processing inevitably interacts with the surrounding 
environmental system (or the thermal reservoir), which induces the quantum world into 
classical world [TKIHEl]. Thus, it is an important subject to analyze the entanglement 
decay induced by the unavoidable interaction with the environment [H [91 [101 [Hi [12] . In 
one-party quantum system, this process is called decoherence [T3 | [T ^ [T5l [T6l \T7\ [T8 | [T9] . 
In this paper, we will analyze the entanglement dynamics of bipartite non-Markovian 
quantum system. As well known, the system can only couple to a few environmental 
degrees of freedom for short times. These will act as memory. In short time scales 
environmental memory effects always appear in experiments [20]. The characteristic 
time scales become comparable with the reservoir correlation time in various cases, 
especially in high-speed communication. Then an exactly analytic description of the 
open quantum system dynamic is needed, such as quantum Brownian motion(QBM) 
[2~T] . a two-level atom interacting with a thermal reservoir with Lorentzian spectral 
density [22], and the devices based on solid state [23] where memory effects are 
typically non- negligible. Due to its fundamental importance in quantum information 
processing and quantum computation, non-Markovian quantum dissipative systems have 
attracted much attention in recent years [3 (2U [251 [261, 123, [28]. Recently, researches 
on quantum coherence and entanglement influenced and degraded by the external 
environment become more and more popular, most of the works contributed to extend 
the open quantum theory beyond the Markovian approximation [29], [301 IS]- m [29] . 
two harmonic oscillators in the quantum domain were studied and their entanglement 
evolution investigated with the influence of thermal environments. In [30], the dynamics 
of bipartite Gaussian states in a non-Markovian noisy channel were analyzed. All in 
all, non-Markovian features of system-reservoirs interaction have made great progress, 
but the theory is far from completion, especially how the non-Markovian environmental 
influence the system and what the difference is between Markovian and non-Markovian 
system evolution are not clear. 

In this paper we will compare the non-Markovian entanglement dynamics with 
the Markovian one [32] in Ohmic reservoir with Lorentz-Drude regularization in the 
following three conditions: u <C lo c , ujq ps uj c and ujq ^> u c , where Uo is the 
characteristic frequency of the quantum system of interest and uj c the cut-off frequency 
of Ohmic reservoir. Thus, u c <C uoq implies that the spectrum of the reservoir does not 
completely overlap with the frequency of the system oscillator and u ^> uo c implies 
the converse case. Another point of the entanglement dynamics is the temperature. 
We characterize our system by low temperature, /c#T = 0.03a;o, medium temperature, 
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ksT = 3co>o, and the high temperature k B T = 300lo . We give stationary analysis of the 
concurrence [9] and find that the dynamics of non-markovian entanglement concurrence 
C at temperature k B T = is fundamentally different from the k B T > 0. We find 
that "entanglement sudden death" (ESD) depends on the initial state when k B T = 0, 
otherwise the concurrence always disappear at finite time when k B T > 0, which means 
that ESD must happen. Maniscalco S et.al studied the separability function S(r) in 
|30j . where the entanglement oscillation appears for twin-beam state in non-Markovian 
channels for high temperature reservoirs. The main result of this paper is that the 
non-Markovian entanglement dynamics is fundamentally different from the Markovian 
one. In the Markovian channel, entanglement decays exponentially and vanishes only 
asymptotically, but in the non-Markovian channel the concurrence C p (t) oscillates, 
especially in the high temperature case. 

The paper is organized as follows. We first introduce the open quantum system and 
the non-Markovian quantum master equation for driven open quantum systems by the 
noise and dissipation kernels. In Sec. Ill we introduce the Wootters' concurrence and 
the initial "X" states. By substituting the initial states into the master equation we get 
the first order coupled differential equations, and give the stationary analysis. In Sec. 
IV, we numerically analyze the Markovian and non-Markovian entanglement dynamics. 
Then an open-loop controller adjusted by the temperature is proposed to control the 
entanglement and prolong the ESD time. Conclusions and prospective views are given 
in Sec. V. 



2. The model 



Our system consists of a pair of two-level atoms (two qubits) equally and resonantly, 
coupled to a single cavity mode, with the same coupling strength. The master 
equation for the reduced density matrix p{t) which describes its dynamics is given by 

ESI EH ED [33] 

dp{t) A(t)+ 7 (t) ^ ro + + + 

= 2 ^ { j pa ? ~ **** 9 ~ P * j } 

+ m_ 2 vij :{2atpa -_ a - a> _ paM] (1) 

where a + = \{p\ + icrz), cr~ = \(o~i — «cr 2 ), with <ti, a 2 the Pauli matrices. The time 
dependent coefficients appearing in the master equation can be written, to the second 
order in the coupling strength, as follows 

A(t) = drk(r) cos(uot), 



lit) = J dTp(r)sm(uj T) 

with 



(2) 



k{r) = 2 J °° du J(ui) coth[ui/2k B T} cos(ur), 
Ii{t) = 2 J °° dcu J(co) sin(o;r), 
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being the noise and the dissipation kernels, respectively. This master equation (1) 
is valid for arbitrary temperature. The coefficient 7(2) gives rise to a time dependent 
damping term, while A(t) the diffusive term. The non-Markovian character is contained 
in the time-dependent coefficients, which contain all the information about the short 
time system- reservoir correlations [7j. In the previous equations J(u) is the spectral 
density characterizing the bath, 

'M = f£^(— ) M 

i 

and the index i labels the different field mode of the reservoir with frequency tUj. Let 
the Ohmic spectral density with a Lorentz-Drude cutoff function, 

7T LO A c + UJ Z 

where u is the frequency of the bath, and uj c is the high-frequency cutoff. 
Then the closed analytic expressions for A{t) and j(t) are [TBI I3T] 

2 

1® = T^l 1 ~ cos ^t) - re~ r ^ sm(u t)}, (6) 

1 + r z 

A(t) = ujq ^{coth(7rr ) — cot(7rr c )e _a;ci [r cos(u t) — sin(u; t)] 

1 + r 1 

+ — cos(^ot) [F(-r c , t) + F(r c , t) - F(ir , t) - F(-ir , t)} 
nr 

1 e~ Ult 
-sin(a;ot)[— — 
7T 2r (l+rg) 

+^F(-r c ,t) - F(r c ,t)]\h (7) 

where r = u; /27r/csT, r c = ou c /2irkBT, r = u c /u , and 

F{x,t) = 2 F^x, 1,1 + x,e- vit ), (8) 

^,1)^^(2,1 + 1,2 + 1,6^). (9) 

v\ = 2nkBT, and 2-^1(0, b, c, z) is the hypergeometric function. Note that, for time t 
large enough, the coefficients A(t) and j(t) can be approximated by their Markovian 
stationary values Am = A(t — > 00) and 7m = j(t — * 00). From Eqs. (6) and (7) we 
have 

„2 



■-Bm(u t)[ — ^ ; — ^r[(r - i)G(-r , t) + (r + i)G(r ,t)] 



im = +++T, (10) 

1 + r z 



and 



j .2 



A M = UJ o l + r 2 coth(7rr ). (11) 

Note that 7(2) has nothing to do with the temperature [33] . In Fig.l we plot the time 
evolution of non-Markovian coefficients A(t) and 7(2) in different channel temperatures. 
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Figure 1. (Color online) Dynamics of non-Markovian coefficients A(t) (blue solid line) 
and "f(t) (red dotted line) at different temperatures: (a)low temperature fcgT = O.Of , 
(b)medium temperature fcgT = 1, and (c)high temperature ks{t) = 100, respectively. 
The other coefficients are chosen as r = 0.1 ujq = 1, and a 2 = 0.01. 



In Fig. 1(a), the temperature is k^T = 0.01. There are two important main points 
embodied in the Figure, the first is that the coefficient 7(t) has dominated the system 
dissipation at low temperature, the other Am = 7a/ i n the long time limit. Fig. 1(b) 
and (c) are the evolution at the medium temperature and high temperature respectively. 
The Figure shows that the larger the temperature, the more important the coefficient 
A(t). 



3. Concurrence and initial states 



In order to describe the entanglement dynamics of the bipartite system, we use the 
Wootters concurrence [9], Elj. For a system described by the density matrix p, the 
concurrence C(p) is 

C(p) =max(0, ^Y 2 ~ y/h - yf>*)> (12) 

where Ai, A 2 , A 3 , and A 4 are the eigenvalues (with Ai the largest one) of the "spin-flipped" 
density operator £, which is defined by 

c = ^^>*«®0. ( 13 ) 

where p* denotes the complex conjugate of p and a y is the Pauli matrix. C ranges in 
magnitude from for a disentanglement state to 1 for a maximally entanglement state. 
The concurrence is related to the entanglement of formation Ef(p) by the following 
relation [31] 

E f (p) = e[C(p)], (14) 
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where 



e[C{p)\ = h[ 



1 + V /1-C 2 (p). 



and 



h(x) = — x\og 2 x — (1 — x) log 2 (l 



•J' 



(15) 



(16) 



Assume that the system is initially an "X" state, which has non-zero elements only 
along the main diagonal and anti-diagonal. The general structure of an "X" density 
matrix is as follows 

Pl4 \ 





P 





P22 
P23 





P23 
P33 



(17) 



/ Pu 


P* 23 P33 

\ Ph p 4 4 / 

Such states are general enough to include states such as the Werner states, the maximally 
entangled mixed states (MEMSs) and the Bell states; and it also arises in a wide variety 
of physical situations [351 ESI EZ]- This particular form of the density matrix allows us 
to analytically express the concurrence as [38] 



Cf = 2mgx{0,K 1 ,K 2 }, 



where 



K 2 



\P23\ 



y/PllpU, 



(19) 



y/p22p33- 

A remarkable aspect of the "X" states is that the time evolution of the master equation 
(1) is maintained during the evolution. Substituting (17) into (1), the non-markovian 
master equation of two-qubits system, we obtain the following first-order coupled 
differential equations: 

p n (t) = -2(A(t) + 7(<))Pii(*) + (A(t) - 7(i))P22 (f) + (A(t) - -f(t))p 33 (t), 
p 22 (t) = (A(t)+ 1 (t)) Pll (t)-2A(t)p 22 (t) + (A(t)- 1 (t))p 44 (t), 
p 33 (t) = (A(t)+y(t))pn(t)-2A(t)p 33 (t) + (A(t)-i(t))pu(t), 
p M (t) = (A(t) + i(t))p 22 (t) + (A(t) + j(t))p 33 (t) - 2(A(t) - *y(t)) P4A (t), 
p 23 {t) = -2A{t)p 23 {t), 
p u (t) = -2A(t)p u (t). 

From Eq. (18) the concurrence C is dependent on the coefficients A(t — > oo) and 
j(t — > oo) in the asymptotic long time limit. Eqs. (10) and (11) give the stationary 
value of 7(t) and A(t), the Markovian limit 



(20) 



7a/ = y(t -> oo) 



and 



A{t 



oo 



1 + r 2 



coth( 



2k B T' 
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7m doesn't depend on temperature, but A M is monotonically increasing with respect 
to temperature T. When T -> 0, A M -> f^. Noting coth(vrr ) ~ 1 + ^ ~ ^f, at 
high temperature 

A £ T = ^bTj^- (21) 

So Am > 7j\/ is noticeable when temperature ksT > 0. From Eqs. (20) we can get the 
stationary solution 



Pii(*-«0 = f^p 33 (t-oo) 
p 22 (t -> oo) = p 33 (t ->■ oo), 
p 33 (t -> oo) = 



AA' 2 ' 



(22) 



and 



p 23 (t -> oo) = 0, 
Pi4(£ - * oo) = 0. 
According to Eqs. (18, 19), 



(23) 



A 2 2 

i^, 2 (t - oo) = - m-7at < (24) 

This means that entanglement must disappear in a finite time period, i.e. the ESD must 
happen. 

When temperature ksT = 0, Am ~ 7a/- From Eqs. (20) we can also get the 



(25) 



stationary solution 






Pu(t - 


-> oo) 


= 


P22{t ~ 


-> oo) 


= 


P33(t - 


-> oo) 


= 


p AA {t - 


-> oo) 


= 1 


and 






P2z{t ~ 


■» oo) 


= 


Pu(t - 


-> oo) 


= 



(26) 

— * cxj j = u. 

From Eqs. (18, 19), 

K lj2 (t oo) = 0. (27) 

This means that entanglement maybe disappear asymptotically, or oscillates, or other 
complex behaviors. In the following, we use the numerical methods to demonstrate the 
concurrence evolution for a special kind of "X" state, the pye state. 
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r=0.1 low temperature, non-Markovian r=1 low temperature, non— Markovian r=10 low temperature, non— Markovian 




Figure 2. (Color online)Time evolution of non-Markovian concurrence as a function 
of parameter "a" in the low temperature reservoirs. 



r=0.1 medium temperature, non-Markovian r=1 medium temperature, non-Markovian r=10 medium temperature, non-Markovian 




Figure 3. (Color online)Time evolution of non-Markovian concurrence as a function 
of parameter "a" in the medium temperature reservoirs. 
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Figure 4. (Color online)Time evolution of K{t) for temperature ksT = O.OOOOOlwcb 
r = 0.1, and initial state pye for the cases a = (black dotted), a = 0.3 (cyan 
dash-dotted line), a = 0.5 (red dash line), a = 0.6 (green dotted-dotted line), a = 0.8 
(magenta asterisk), and a = 1.0 (blue solid line). 



4. Non-Markovian vs. Markovian entanglement dynamics 



In this section, we use the formalism of the preceding section to determine the 
disentanglement. As an example, let us consider an important class of mixed states 
with a single parameter a like the following [23, [39J HO] 



/ 



Pye 



V 




1 
1 








1 — a 



\ 



J 



Apparently, the concurrence of pye is C p (t) = max{0, K(t)}, and K(t) = 

|[1 — \J a(l — a)}. In our simulations, 



^Pu(t)pu(t). Initially, C(p(0)) 



(28) 

P23(t)\ ~ 
UJ = 1 

are chosen as the norm unit, and we regard the temperature as a key factor in 
disentanglement process, for high temperature /c#T = 300a;o, intermediate temperature 
k B T = 3u , and low temperature k B T = 0.03a;o, respectively. Another reservoir 
parameter playing a key role in the dynamics of the system is the ratio r = uj c /ujq 
between the reservoir cutoff frequency uj c and the system oscillator frequency uoq. As we 
will see in this section, by varying these two parameters k B T and r = uj c /uj , the time 
evolution of the open system varies prominently from Markovian to non-Markovian. 

In Fig. 2, the time evolutions of the non-Markovian concurrence for various values 
of the parameter a in low temperature is plotted. From Fig. 2 we can see that the 
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r=0.1 high temperature, non-Markovian r=1 high temperature, non-Markovian r=10 high temperature, non-Markovian 




r=0.1 high temperature, Markovian evolution r= i high temperature, Markovian evolution r=10 high temperature, Markovian evolution 




Figure 5. (Color online) Comparing the non-Markovian entanglement dynamics with 
the Markovian one by the time evolution of concurrence as a function of parameter 
"a" in high temperature reservoirs, at r = 0.1, r = 1, r = 10 respectively. 



entanglement dynamic relies on the different values of r = u c /uj . If the spectrum of 
the reservoir does not completely overlap with the frequency of the system oscillator, 
r < 1, we can see from Fig. 2 that the ESD time is considerably long. As increases 
the ratio r, the ESD time becomes shorter and shorter. With different initial state 
we can see that the concurrence varies prominently. When the initial state a = 0, 
the non-Markovian entanglement decay slowly, as increasing a, the entanglement decay 
intensely, which means that we can prepare certain initial entanglement states and use 
this fact to control the system environment in order to prolong the entanglement time. 

Fig. 3 is the medium temperature case. Like Fig.2, under different systems, different 
entanglement initial states, corresponding to different values of a, and different r, some 
decay faster, some slower. But there are some fundamental difference between Fig.2 and 
Fig. 3. In Section III, we get the concurrence in the long time limit, and we affirmed that 
when temperature ksT = 0, the dynamics of non- markovian entanglement concurrence 
C is fundamentally different from the case of ksT > 0. As we can see from Fig. 3, 
for "py^" states, as soon as the temperature larger than zero, the concurrence always 
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r=0.1 a=0, non— Markovian entanglement evolution r=0.1 a=0, Markovian entanglement evolution 




k B T o> t k B T co Q t 



Figure 6. (Color online) Comparing the non-Markovian entanglement dynamics with 
the Markovian one by the time evolution of C p (t) as a function of "ksT" for initial 
state a = and r = 0.1. 

disappear at finite time and there were no long-lived entanglement for any value of a, 
which means that ESD must happen. The theoretical proof is K(t — * oo) < 0. But 
when = 0, the stationary value of Kit — > oo) equals zero. So, whether or not 
and when the ESD will happen are not sure in k^T = 0. In Fig.4 we give a numerical 
analysis of entanglement dynamic with different initial states and find that there exists 
a £ G (0, 1), for almost all values a > £, the concurrence is completely vanished at a 
finite time, which is the effect of ESD. However, for < a < £, the entanglement of 
this state decays exponentially. But when t — > oo, for all initial state, i.e. a G [0, 1] the 
concurrence will tend to be 0. 

Fig. 5 is the high temperature case. One of the remarkable phenomenon in this 
figure is that the ESD time is short. In typical experimental conditions, quantum 
dots are subjected to an external magnetic field B ~ 1 — 10T jl6], the ESD time 
tESD ~ (3 x lCT 1 — 3)//cbT. Another obvious phenomenon is in high temperature the 
Markovian quantum system decays exponentially and vanish only asymptotically, but in 
the non-Markovian system the concurrence C p (t) oscillates, which is evidently different 
from the Markovian. In this case the non-Markovian property becomes evidently. This 
oscillatory phenomenon is induced by the memory effects, which allows the two qubit 
entanglement to reappear after a dark period of time. This phenomenon of revival of 
entanglement after finite periods of "entanglement death" appears to be linked to the 
environment single qubit non-Markovian dynamics, in particular, the A(£) < at some 
times in some environment [31]. The physical conditions examined here are, moreover, 
more similar to those typically considered in quantum computation, where qubits are 
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taken to be independent and where qubits interact with non-Markovian environments 
typical of solid state microdevices [4"T] . 

As we indicated above, temperature is one of the key factor in the entanglement 
dynamic. Figs. 2, 3, 4, 5 are plotted in the chosen temperature, while in Fig. 6 
k B T ranges from to 100. In Fig. 6 the concurrence vs "temperature k B T n vs u^t in 
r = 0.1, and the initial state is the u Xye" state with a = 0. From Fig. 6 we can compare 
the non-Markovian entanglement dynamics with the Markovian one clearly. The left 
is the non-Markovian one from which we can see the oscillation of the concurrence. 
Moreover, at the temperature the non-Markovian effect is faint, as the temperature 
rises, the non-Markovian becomes more and more obvious, while the Markovian one 
decays exponentially. This phenomenon embodies the non-Markovian effect, which 
is evidently different from the Markovian property. Maniscalco S et. al studied the 
separability function S(t) in [30], where entanglement oscillation appears for twin-beam 
state in non-Markovian channels in high temperature reservoirs. Both of them have the 
same phenomenon. Ref. [31] gave a distribution curve when A(r, t) — j(r,t) > and 
A(r, t) — j(r,t) < 0. We convince that due to the non-Markovian memory effect, 
particularly A(t) < in Eqs (20), the entanglement concurrence oscillates. With 
A(t) — 7(t) > the concurrence descended whilst A(t) — -f(t) < the concurrence 
ascended, which guide us to adjust the temperature to control the entanglement 
evolution. In order to show this and motivate the related research we design the open 
loop controller 

k B T = e- alA{t) -^ t)l k B T (29) 

where a is the modulation, and k B To is the initial temperature. In Fig. 7, we plot the 
controlled entanglement evolution, where the initial temperature is chose as k B T = 30, 
which oscillates and ESD occurs at t ~ 19. According to Fig. 1, 7(t) can be neglected. 
For different modulation a, different controlled entanglement evolution is plotted, and 
the ESD time can be prolonged for considerable long time. 

5. Conclusions 

In this paper we have presented a procedure that allows to obtain the dynamic of a 
system consisting of two identical independent qubits, each of them locally interacting 
with a bosonic reservoir. A non-Markovian master equation between two qubit systems 
in the same environment was obtained. We characterize our entanglement by the 
temperature and the ratio r between ujq the characteristic frequency of the quantum 
system of interest, and u c the cut-off frequency of Ohmic reservoir. For a broad class 
of initially entangled states, "X" states, by useing Wootters' concurrence, we analyze 
the long time limit phenomenon of the entanglement dynamic. We find that the 
dynamic of non- markovian entanglement concurrence C p (t) at temperature k B T = 
is fundamentally different from k B T > 0. When k B T = 0, from our numerical 
analysis, we find that "entanglement sudden death" (ESD) occurs depending on the 
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Figure 7. (Color online) Controlled entanglement evolution with different modulation 
a = 3 (green dotted line), a = 2 (red dashed line), a = 1 (blue dashed-dotted line), 
and initially evolution (black solid line), respectively. 

initial state, but if ksT > the concurrence always disappear at finite time, which 
means that ESD must happen. In the ksT = case, we find that there exist a 
£ G (0,1), for all values a > £, the concurrence is completely vanished in a finite 
time, which is the effect of ESD. However, for < a < £, the entanglement of this 
state decays exponentially. But when t — > oo, for all initial state, i.e. a e [0, 1] 
the concurrence will tend to be 0. From our numerical analysis we also find that the 
entanglement dynamic relies on the different values of r = u c /ujq. If r <C 1, the ESD 
time is considerably long. As increases the ratio r, the ESD time becomes shorter 
and shorter. Moreover, when the initial state a = 0, the non-Markovian entanglement 
decays slowly, as increases a, the entanglement decays intensely. Most of all, we have 
shown that the non-Markovian dynamics of entanglement, described by concurrence, 
presents oscillation even revivals after entanglement disappearance, typically for high 
temperature non-Markovian system. At last, we design the open loop controller which 
adjust the temperature to control the entanglement and prolong the ESD time. 
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